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Abstract 

The quantum H<± integrable system is a 4D system with rational potential related to the non- 
crystallographic root system H4 with 600-cell symmetry. It is shown that the gauge-rotated H4 
Hamiltonian as well as one of the integrals, when written in terms of the invariants of the Coxeter 
group H^, is in algebraic form: it has polynomial coefficients in front of derivatives. Any eigen- 
functions is a polynomial multiplied by ground-state function (factorization property). Spectra 
corresponds to one of the anisotropic harmonic oscillator. The Hamiltonian has infinitely-many 
finite-dimensional invariant subspaces in polynomials, they form the infinite flag with the charac- 
teristic vector a = (1,5,8,12). 
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I. THE HAMILTONIAN 



k=l L 



The quantum H± system is a four-dimensional system related to the non-crystallographic 
root system if 4 The Hamiltonian of this model is invariant with respect to the H 4 
Coxeter group, which is the full symmetry group of the 600-cell polytope. The H4 Coxeter 
group is discrete subgroup of 0(4) and its dimension is 14400. In Cartesian coordinates the 
if 4 rational Hamiltonian has the form (see 

9 ?) T 2 U X k x 2 [ + (_1)M2 X2 _|_ (-1)W X + (-1)/M X4 ]2 

K K - > M2, 3, 4=0,1 ^ 

+ v V ^ 

where = {1,2,3,4} and its even permutations. Here g = v{y — 1) > —1/4 is 

the coupling constant, ip± = (1 ± v / 5)/2 the golden ratio and its algebraic conjugate. The 
configuration space is the subspace of R 4 where the condition (a-x) > holds for any positive 
root a of if 4. It is an analogue of the principal Weyl chamber in the case of crystallographic 
root systems. 

The ground state eigenfunction and its eigenvalue are 

ttoOO = A? A£A* ex P ( " I E x l) > E o = Ml + 30l/ ) ' ( 2 ) 



where 



i 



Ai = ] Jx fcl (3) 



fc=i 



a 2 = n [x 1 + (-irx 2 + (-irx 3 + (-1)^x4], (4) 

M2, 3, 4=0,1 

a 3 = n n ^ + (-1)^+^ + (-ir^it + ■ x,] . (5) 

{ij',fe,0 Mi, 2=0,1 

The ground state eigenfunction ([2]) does not vanish in the configuration space. 

The main object of our study is the gauge-rotated Hamiltonian with the ground state 
eigenfunction (j2J) taken factor, 

V = ^(^o)- 1 ^-^)^) , (6) 

where E is the ground state energy given by (J2J). The gauge rotated operator <^ is the 
second-order differential operator without free term. By construction its lowest eigenfunction 
is a constant and the lowest eigenvalue is equal to zero. 



Now let us define new variables in ([6]). The if 4 root space is characterized by four 
fundamental weights w c , c = 1,2,3,4 (see e.g. ^j). Taking action of all group elements on 
fundamental weight u c we generate orbit Q(w c ) of a certain length (length = ^elements of 
the orbit). The results are summarized as 





weight 


orbit length 


Wi 


= (0, 0, 0, 2^+) 


120 




= (1, ft, 0, ft) 


600 


w 3 = 


(0, ip + , 1, ft - 1) 


720 


= 


= (0, 2tp + , 0, 2ft) 


1200 



Now let us find ^-invariants. In order to do it we choose for simplicity the shortest orbit 
Q(wi) and make averaging, 

€ } (x)= £ (w-x)\ (7) 

where a = 2,12,20,30 are the degrees of the H4 group. It is worth noting that these 
invariants are defined ambiguously, up to a non-linear combination of the invariants of the 
lower degrees 
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where {A} are parameters. Canonical invariant basis for the H4 was found only recently by 
Mehta j4] (see also [ij]). Now we use if 4 invariants as new coordinates in 

Now we can make a change of variables in the gauge- rotated Hamiltonian ([6]): 

( T T T T \ _v (d n ) 

^Xi, J/2; -^3) x 4) ^ 2 ) t 12 ) fc 20 ) ''SO J • 

The first observation is that the transformed Hamiltonian hn 4 (t) ([6]) takes on an algebraic 
form for any value of the parameters {^4} in variables t's ([8]). The second observation is 
that for any value of the parameters {A} the operator hn 4 (t) has infinitely- many finite- 
dimensional invariant subspaces in polynomials which form infinite flag. Our goal is to find 
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the parameters for which hn 4 (t) preserves a minimal flag (for a discussion see e.g. [3]). After 
some analysis such a set of parameters is found 



43510 „ 41701 „ 17583778485 

Ai = — 1, A 2 = , A 3 = , At = 

1809 1809 146142376 

313009515 „ 22081114965 „ 798259915667 

As = ~ irooo ^ , A 6 = —— — — — , A 7 



(9) 



15383408 7691704 ' ' 292284752 

Hereafter we call the t-variables for such values of parameters as r-variables. 

In order to write down explicit expressions for variables r it is useful to exploit the 

notation for multivariate polynomials introduced by Iwasaki et al For given partition 

A = (Ai, A 2 , A3, A 4 ) with Ai > A 2 > A 3 > A 4 > 0, let denote as M x the associated monomial 

symmetric polynomial of the variables (x 2 , x\, x\, a; 4 ), 



2fi 2 2fi 3 2(j, 4 

X o Xo Ju a 



where the sum is taken over all permutations (/ii, /12, ^3, of (Ai, A2, A3, A4). If A consists 
of mutually distinct numbers p\ > ■ • • > p m with pj appearing kj times in A, then we denote 
the polynomial as 

Mx = W I • ■ ■ \p k -] ■ 

Let introduce also A 4 as the fundamental alternating polynomial of (xf, x|, x\, x 4 ), 



a 1= n 

l<i<j<4 

Then, the r-variables (7)-(9)are written in these notations as 



n = [l|0 d ] = x( + x\ + x\ + x\ , 



r 2 = 14[3 2 |0 2 ] -6[4|2|0 2 ] +2[5|1|0 2 ] - 270[2 2 |1 2 ] + 30[2 3 |0] 
-12[4|1 2 |0] + 348 [3|1 3 ] + 9[3|2|1|0] + 33v / 5A 4 , 

r 3 = 2[8|2|0 2 ] +4[8|1 2 |0] - 10[7|3|0 2 ] -45[7|2|1|0] + 60[7|1 3 ] 
+22[6|4|0 2 ] + 157[6|3|1|0] + 270[6|2 2 |0] - 150[6|2|1 2 ] 
-22[5 2 |0 2 ] - 131[5|4|1|0] - 733[5|3|2|0] - 2156[5|3|1 2 ] 
+4050[5|2 2 |1] + 1320[4 2 |2|0] + 4650[4 2 |1 2 ] + 6[4|3 2 |0] 
-2175[4|3|2|1] - 19050[4|2 3 ] + 10800[3 2 |2 2 ] + 3336[3 3 |1] 
+3v / 5A 4 { 5[4|0 3 ] - 18[3|1|0 2 ] + 49[2 2 |0 2 ] + 3[2|1 2 |0] 
+1146[1 4 ] } , 



r 4 = 65742[15|0 3 ] - 504[13|2|0 2 ] + 830[13|1 2 |0] + 61690[12|3|0 2 ] 
-5130[12|2|1|0] - 9495[12|1 3 ] + 18795[11|4|0 2 ] 
+28560[11|3|1|0] - 43500[11|2 2 |0] - 53070[11|2|1 2 ] 
-156330[10|5|0 2 ] + 59130[10|4|1|0] + 26415 [10 1 3 1 2 1 0] 
+405255[10|3|1 2 ] + 1350 [10 1 2 2 1 1] + 19710[9|6|0 2 ] 
-20[9|4|2|0] - 8663355[9|4|1 2 ] - 120[9|3 2 |0] + 450[9|3|2|1] 
-962715[9|2 3 ] + 13860 [8 1 7| 2 ] - 94530 [8 1 6 1 1 1 0] 
-353160[8|5|2|0] - 1452060[8|5|1 2 ] + 5557050 [8 14| 3 1 0] 
+590580[8|4|2|1] - 198270 [8|3 2 |1] + 389250 [7 2 |1|0] 

+2897820[7|6|2|0] - 5227920 [7|6| l 2 ] + 1134540[7|5|3|0] (10) 
-4041270[7|5|2|1] - 591330[7|4 2 |0] + 23417850[7|4|3|1] 
-22770[7|4|2 2 ] - 23528790[7|3 2 |2] + 29647380 [6 2 |3|0] 
+36597510[6 2 |2|1] - 1649925 [6 1 5 14 1 0] + 150[6|5|3|1] 
+40935[6|5|2 2 ] - 510[6|4 2 |1] - 60[6|4|3|2] + 242505[6|3 3 ] 
+270060[5 3 |0] - 528270[5 2 |4|1] - 36255[5|4 2 |2] + 825[5|4|3 2 ] 
+707085[4 3 |3] +45v / 5A 4 { -27040 [9|0 3 ] -5[8|1|0 2 ] 
-1914[7|1 2 |0] + 23[6|3|0 2 ] + 91[6|2|1|0] - 44[6|1 3 ] 
-352[5|4|0 2 ] + 8[5|3|1|0] + 1085[5|2 2 |0] + 6875[5|2|1 2 ] 
-5168[4 2 |1|0] - 934 [4|3|2|0] - 568[4|3|1 2 ] + 1773 [4|2 2 |1] 
+20911[3 3 |0] + 15915[3 2 |2|1] + 573[3|2 3 ] } . 



Thus, the variables 7"i,2,3,4 are homogeneous polynomials in x 2 's of the degrees 1,3,10,15 , 
respectively. 

Finally, the gauge-rotated Hamiltonian ([6]) in the r-coordinates (fTUl) written as 



takes amazingly simple form with the coefficient functions 




(11) 



A 



u = 4 ri , A 



24 r 2 , 



40 r 3 , A 



60 r 4 , 
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A 22 = 88 nr 3 + 8 t\t 2 , 

A 23 = -4 t\t\ + 24 rfrs - 8 r 4 , 

A 24 = 10 r^rl + 60 rfr 2 r 3 + 40 t\t± - 600 r 3 2 , 

38 8 

-433 = -y nrf + 28 r^Ta - - r 1 V 4 , 

A34 = 210 Tfrfra + 60 r 1 3 r 2 r 4 - 180 t 4 t 2 + 30 r 2 4 , 

A 44 = -2175 rir 2 3 r 3 - 450t 1 2 t 2 2 t 4 - 1350 t\t 2 tI - 600 r 1 4 r 3 r 4 , (12) 



B x = 8(1 + 30z/) - 4un , 

B 2 = 12(1 + Wu) rf - 2Aujt 2 , 

B 3 = 20(1 + 6u) t\t 2 - 40wr 3 , 

B A = 15(1 - 30v) t\tI - 450(1 + 2u) r 4 r 3 - 60u;r 4 . 

It can be easily checked that the operator ( fTTT) is triangular with respect to action on 
monomials rf 1 r 2 2 rf 3 r 4 4 . One can find the spectrum of ( TTTT) h^ip = —2e(p explicitly 

em,n 2 ,n3,n4 = 2w(rii + 6n 2 + 10n 3 + 15n 4 ) , (13) 

where rij = 0, 1, 2, . . .. Degeneracy of the spectrum is related to the number of solutions of 
the equation n\ + 6n 2 + 10n 3 + 15n 4 = n for n — 0, 1, 2 ... in non-negative numbers ni )2i3)4 . 
The spectrum e does not depend on the coupling constant g and it is equidistant. It coincides 
to the spectrum of 4D anisotropic harmonic oscillator with frequencies (2a;, 12a;, 20a;, 30a;). 
The energies of the original rational if 4 Hamiltonian ([T|) are E = E + e. It is worth noting 
that the Hamiltonian ( TTTT) has infinite family of eigenfunctions ni ,0,0.0 depending on single 
variable T\. They are given by the Laguerre polynomials and the eigenvalues are linear in 
quantum number (cf. (fT3l) ) 

^,0,0,0^1) = L^ +mu \un) , e ni , ,o,o = 2um , m = 0, 1, 2, . . . . (14) 

The boundary of the configuration space of the rational if 4 model ([I]) in the r variables is 
determined by the zeros of the ground state eigenfunction, hence, by pre-exponential factor 
in (T5]). It is the algebraic surface of degree 120 in Cartesian coordinates being a product of 
monomials. In r-coordinates (flOl) it can be written as 
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64 t 15 t 4 3 + 1440 t1 4 t 2 t z tI + 10800 t^t^t^ + 27000 t 12 t 3 t 3 - 240 t 12 t 3 t 4 2 

- 3600 r 1 11 r 2 4 r 3 r4 - 13500 t 10 t 2 5 t 3 2 + 34992 t^t\ - 1440 r^r 2 + 300 r 1 9 r 2 6 r 4 

- 2160 rjVsrf r 4 - 144(hfr 2 7^ + 2250 tJt 3 - 22680 7^ - 28080 rf^r 2 

- 203760 r 1 7 T 2 3 T 3 2 T4 - 125 rfrf - 493020 r^rf + 3600 ?f t|t| + 57780 r 1 5 r 2 5 r 3 r4 

- 8640 r 1 5 r 3 4 r4 + 4320 r^rl + 221310 t 4 t 2 6 t 3 2 - 648000 r^r 2 r| + 116640 r^r^r 2 

- 4680 r 1 3 r 2 7 r 4 + 712800 r 1 3 r 2 2 r 3 3 r4 + 6480 r 3 r 2 2 r 3 - 35640 t 2 t 8 t 3 + 2052000 rf^r 4 

+ 62640 Tfrfrgr 2 + 259200 TiT-frf r 4 + 1944 r 2 10 + 129600 r|r 3 3 + 2592 r 2 5 r 4 2 

(15) 

+ 2160000 r| - 86400 r 3 r 2 + 864 r 4 4 = , 
which is the algebraic surface of degree 18 being given by a polynomial of degree 15 in n, 
of the degree 10 in r 2 , of the degree 6 in r 3 and of the degree 4 in t±. It is worth mentioning 
that l.h.s. of (fT5|) is proportional to the square of Jacobian, </ 2 (§^)- 

The Hamiltonian h-^ir) has infinitely- many finite-dimensional invariant subspaces 

P (i,5,8,i2) = ^ T P W P2^ 3 ^ \o< Pl + 5p 2 + 8p 3 + 12p 4 < n) , n = 0, 1, 2, . . . , (16) 
which form the (minimal) infinite flag. Its characteristic vector is 

a min = (1,5,8,12) . (17) 

It is worth noting that each particular space •p^ 1 ' 5,8 ' 12 ) ([16]) as well as the whole flag are 
invariant with respect to a weighted projective transformation 

Ti -»■ ri + a , 

r 2 ->■ r 2 + &i rf + 6 2 r x 4 + 6 3 rf + 6 4 r x 2 + b 5 n + 6 6 , 

T3 -> r 3 + ci r 3 r 2 + c 2 r 2 r 2 + c 3 rir 2 + c 4 r 2 + c 5 r 8 + c 6 r x 7 

+c 7 rf + c 8 rf + c 9 r x 4 + cio rf + c n rf + c i2 n + c 13 , (18) 
r 4 ->■ r 4 + di r 4 r 3 + d 2 r 3 r 3 + d 3 r 2 r 3 + d 4 TiT 3 + d 5 t 3 
+d 6 t\t 2 + d 7 t\t 2 + d 8 r 2 + d 9 t\t 2 + d 10 r 3 r 2 
+dn r 2 r 2 + di2 nr 2 + d 13 r 2 + d u r 12 + d 15 r 11 
+dx e r 10 + di7 if + d 18 r 8 + di 9 + d 20 if + d 21 r\ 
+d 22 r 4 + d 23 r 3 + d 24 r 2 + d 25 n + d 26 , 

7 



where {a, b, c, d} are parameters. It manifests a hidden invariance of the Hamiltonian ([I]) 
preserving its algebraic form. A meaning of this invariance is unclear. 



II. INTEGRAL 

The Hamiltonian (JT]) can be written in hyperspherical coordinates 

x\ = r sin ij) sin 9 cos <fi , 
X2 = r sin ip sin 9 sin <fi , 
£3 = r sin?/; cos 9 , 
£4 = rcos^ , 



(19) 



where it takes the form 



„^_I AW + I W V + M). (20) 



Here A^ 4 ^ is the AD Laplacian and the angular function 



a ^ 2v(u-l) 2v{v 



2v(v-l) 2v(v-l) 



2 



+ 
+ 



2v(u-l) 2v(v-l) 



2i/(i/ -1) 2i/(i/ — 1) 



2 



2v(v-\) 2v(v-l) 
+ 7 " " rrr + 



{s^S B C^ + (f+C^ + (ps^sos^) 2 {s^sqc^ - ip_c^, + ip_s 1 ps e s ( p) 

2u(v-l) 2v{v-\ 
+ 1 : ^ + 



2 



+ 



(s^SeCf, + <p + ty - (p-SipSgScf,) 2 (s^SeCf, - (p + C^ - (p.S^SgS^) 2 

2v(v-l) 2v(v-l) 



(s^SgStt, + (f+S^SgCcf, + (fi-C^) 2 (s^SeSj, - (f+S^SeC^ + (f-Cip) 2 

2v(v-l) 2u(u-i) 



+ 



2v{y-\) 2v{v-l) 



2 



2v{v-l) ^ 
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2v{v-l) 2v(u-l) 

(s^sgs^ + ip + c^ + ip-s^ce) 2 (s^sgs^ - (p + c^ + cp-S^cg) 2 
2v(y-l) 2v(v-l) 



(s^sgs^ + v9 + c^ - (p-s^c e ) 2 (s^sgs^ - ip+c^ - ips^cg) 2 
2v(v-l) 2v(v-l) 



(s^cg + ip+SipSec^ + (f-S^sgs^) 2 (s^,cg + (p-S^sgc^ + (ps^sgs^) 
2v{v-l) 2v{v-V) 



2 



(s^pce + (p+s^sgc^ - tp-S^soStp) 2 (s^cg + (ps^sgc^ - (p-s^ses^) 
2v(v-l) 2v(v-l) 



2 



(S^Q? + + (pS^SgC^) 2 (s^Cg - (p + C^ + tpSipSoC^) 2 

2v(v-\) 2v{v-\) 

(s^cg + ip+ty - (f-s^sgc^) 2 (s^cg - (f+c^p - (p^s^sgc^) 2 

2v{v-l) 2v{y-\) 

(s^cg + cp+s^sgs^ + (p-C^) 2 (s^.cg + (ps^sgs^ + (p-C^) 

2v{v-l) 2v{v-\) 



2 



+ 



(s^Cg + If+S^SgS^ - If-C^p) 2 (s^Cg + (pS^SgS^ - y?_C</,) 2 

2v{y-\) 2v{v-l) 

(C^j + If+S^SgCtj, + (fS^Cg) 2 (c^ + (f^S^SgC^ + If-S^Cg) 



2 



I 2v(y-\) | 2v{y-\) 

(c^ + V9 + S^SeC^ - (f-S^Cg) 2 (ty + ifS^SgC^ - ifS^Cg) 2 

2v{y-\) 2v{v-\) 



2v{y-\) 2v{v-V) 



+ 



(ty + (f+S^SgS^ - tpS^SgC^) 2 (c^, - if+S^SgS^ - (f^S^SgC^) 2 

2v{v-V) 2v{v-V) 

(c^ + (fi+S^Cg + ip^S^SgS^) 2 (ty - (fi+S^Cg + ip^S^SgS^) 



2 



2v{y-\) 2v{y-\) 



(C^ + if+S^Cg - (f-SipSgS^) 2 (ty - (p+S^Cg - tpS^SgS^) 

2v{v-l) 2v{y-\) 



2 



+ 
+ 



(s^SgC^ + S^SgS^ + S^Cg + C^) 2 (s^SgC^ - S^SgS^ + S^Cg + C^)'' 

2v{v-l) 2v{y-\) 



(s^SgC^ + S^SgS^ - S^pCg + C^) 2 (s^SgC^ + S^SgS^ + S^Cg - C^)' 

2v{y-\) 2v{v-\) 



(s^SgC^ - S^pSgStp - S^pCg + C^) 2 (s^SgC^ - S^SgS^ + S^Cg - C^,) 5 

2v(v-l) 2v{v-\) 



(s^SgC^ + S^SgScf, - S^pCg - C^) 2 (s^SgC^ - S^SgS^ - S^,Cg - C^) 2 
u{u + \) u(u + l) 1/(1/ + 1) u{u + l) 

' o e 2„2„2 ' n 2 e 2 e 2 o„2„2 9 2 

Here, for the sake of simplicity we denoted = cosip, = sin-^, eg = COS6 1 , = sin^ and 
= cos</>, = sin0. It is seen immediately, that the Schroedinger equation (1201) admits a 



9 



separation of radial variable r: any solution can be written in factorized form 



*(r,iP,9,(f>)=R(r)Q(i>,6,<f>) • 



Functions R and Q are the solutions of the equations 



(22) 



1 d 







, 2 2, 7 

+ -co r + — 



R(r) = ER(r) 



2r 3 Or \ dr J 2 

j-Q(^,M)=7 , 

respectively, while 7 is the constant of separation. The operator J 7 has the form 



(23) 
(24) 



2 sin ip 



d ( . 2 , d \ 1 d ( . . d 
sm ip-^~ I H — : — - — I sine/- 



dtp 



dip 



sin 9 89 



1 «9 2 



<9fl J sin 2 <90 2 



+ (25) 



It can be immediately checked that the Hamiltonian "H# 4 and J 7 commute, 



[U H „F\ = . 



(26) 



Hence, J 7 is an integral of motion. Thus, it has common eigenfunctions with the Hamiltonian 

Let us make a gauge rotation of the operator T (1231) with the ground state function \l/ 
as a gauge factor, 



/ = (^ ) _1 (^ - 7o)^o , 7o = 60i/(l + 30i/) 



(27) 



where 70 is the lowest eigenvalue of J 7 . Then make a change of variables to the r- variables 
(fT0|) . The operator / gets an algebraic form, 

4 o9 4 



/ Yl 7 ''''' ;>-/)--. ' . ' / '''' 7 7 



3=1 



(28) 



where 



^11 = 0, F 12 = 
^13 = 0, F 14 = 



-P22 

^23 
F24 



-4t*t 2 - Urfr 3 + 72r 2 , 

-12t?t 3 + 2r 1 4 r 2 2 + 4nr 4 + 120r 2 r 3 , 

-207fr 4 - 30?f r 2 T 3 - 5rfr 2 3 + 300tit 3 2 + 180r 2 r 4 
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^33 = \t\t A - UrSrs + Y T l T 2 + 200T i > 

F 34 = 90t^ - 30t 1 4 t 2 T4 - 105^^3 - 15rir 2 4 + 300r 3 r 4 , 

2175 

F 44 = 300r 1 5 r 3 r 4 + 675r 4 r 2 r 3 2 + 225t 1 3 t 2 2 t 4 + —r^f r 3 + 450r 4 2 , 

Gi = 0, 

G 2 = -6(1 + lO^ 6 + 12(7 + 60z/)r 2 , 

G 3 = -10(1 + 6z/)r 4 r 2 + 20(11 + 60z/)r 3 , 

G 4 = 225(1 + 2z/)r 1 5 r 3 - y (1 - 30^)^x1 + 40(12 + 45i/)r 4 . (29) 

It is worth noting that in the operator / the variable T\ appears as a parameter. It implies 
that any eigenfunction of the Hamiltonian hn 4 , which depends on n only, is an eigenfunction 
of the integral / with zero eigenvalue. 

It can be also shown that the operator / has infinitely many finite-dimensional invariant 
subspaces in polynomials 

£,(1,6,10,15) = ^i T |2 T |3 T |»4 \o< Pi + q P2 + io p3 + 15 P4 < n) , n = 0, 1, 2, . . . , (30) 

which form a flag with characteristic vector (1,6,10,15). The spectrum of the integral 
J 7 ^/ = T^> can be found in a closed form, 

TototoM = loMMM + 7o = 

72A; 2 + 200fc 2 + 450A; 2 + 120A; 2 A; 3 + 180A; 2 fc4 + 300^4 + 2(l+60z/)(6A; 2 + 10A; 3 + 15fc4)+7o , (31) 

where k 2 , k 3 , fc 4 = 0, 1, 2, . . . and 70 is given by ( )27|) . 

It can be shown that the Hamiltonian /i# 4 has a certain degeneracy - it preserves two 
different flags: one with (minimal) characteristic vector (1,5,8,12) and another one with 
characteristic vector (1,6,10,15). The fact that the operator hn A with coefficients ( fl2l) com- 
mutes with / given by (T28]) implies that common eigenfunctions of the operators /i# 4 and / 
are elements of the flag of spaces , p( 1 ' 6 ' 10 ' 15 ). 

Let us denote <j) n ^ the eigenfunctions of Hh 4 which are elements of the invariant space 
p^i,5,8,i2) their respectful eigenvalues e n ^. The index i numerates these eigenfunctions 
for given n starting from 0. It is evident that an eigenfunction with n < 5 depends on t\ only 
and its 7 is equal to zero. The eigenfunctions with 4 < n < 8 depend on ti j2 only while the 
dependence on the r 12 3 occurs for the eigenfunctions with 7 < n < 12. The eigenfunctions 
with n > 12 depend on all four variables t 1j2)3i4 . 

11 



The function <fi n ^ is related to the eigenfunction of the Hamiltonian "H# 4 ([I]) (and the inte- 
gral J 7 ) through = \l/o0n.i- Thus, the eigenfunctions {0} are orthogonal with the weight 
factor |\l/o| 2 - As an illustration let us give explicit expressions for several eigenfunctions 4>n,i 
and their respectful eigenvalues, 

• n — 

0o,o = 1 , eo,o = , 

• n = 1 

01,0 = ojti - 2(1 + 30z/) , ei j0 = 2w , 

• n = 2 

(j) 2 = u 2 T f - 6w(l + 20z/)ri + 6(1 + 20z/)(l + 30z/) , e 2 , = 4w , 

• n — 3 

3O = w 3 r 1 3 -12w 2 (l + 15z/)r 1 2 + 36w(l + 15z/)(l+20z/)r 1 -24(l + 15z/)(l + 20z/)(l+30z/) , 

£3,0 = 6w , 

• n = 4 

04,o = ^i 1+60iy) (wri) , e 4i0 = 8w , 

• n — 5 

05,o = -4 1+60 ^(wri) , e 5j0 = lOw , 
5il = w 6 r 2 - 3(1 + lOz/VS 5 + 45c/(l + lOz/) 2 ^ 4 - 300w 3 (l + \2v){l + 10i/) 2 r 1 4 + 
900w 2 (l + 15i/)(l + \2v)(l + 10z/) 2 r 2 - 1080w(l + 20z/)(l + 15z/)(l + 12z/)(l + 10i/) 2 Ti+ 
360(1 + 30z/)(l + 20z/)(l + 15i/)(l + 12i/)(l + 10z/) 2 , 
e 5j i = 12w . 

Let us denote ni j the common eigenfunctions of /i# 4 and / which are elements of the 
invariant space p^ 1,6 ' 10 ' 15 ) an d their respectful eigenvalues e n ,iiln,i- The index i numerates 
these eigenfunctions for given n starting from 0. It is evident that an eigenfunction with 
n < 6 depends on t\ only and its 7 is equal to zero. The eigenfunctions with 6 < n < 10 
can depend on and the dependence on the 71,2,3 occurs for the eigenfunctions with 
10 < n < 15. It is evident that all eigenstates (n, i) at fixed n and different i are degenerate: 
their eigenvalues are equal to 2un. We give some eigenstates from p^ 1,6 ' 10 ' 15 ) explicitly, 
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• n = 



0o,o = 1 , eo,o = , 70,0 = 



n = 1 

0i,o = wn - 2(1 + 30v) , ei )0 = 2u , 7 lj0 = 

n = 2 

02,o = u 2 t?-6u(1 + 20v)n + 6(1 + 20u) (1 + 30v) , e 2>0 = 4w , 72i0 = 
n = 3 

3O = w 3 rj 3 -12u; 2 (l + 15z/)r 1 2 +36u;(l + 15z/)(l + 20z/)ri-24(l + 15z/)(l+20z/)(l + 30z/) 

e 3i0 = Qu , 73,0 = 



n = 4 



n = 5 



n = 6 



4,0 = L^ 1+60!/) (wri) , e 4fi = 8u , 7 4> o = 



^5,0 = L£ +e0v '(urri) , e 5:0 = 10w , 75,0 = 



6>o = L^ 1+60l/) (wri) , e 6>0 = 12a; , 7 6 ,o = 
4(7 + oUz/J 

It is worth noting that 05^ = cj 6 06,i + ^06,o where A is a parameter. Eigenfunctions 

0n,O = 0n,O at U < 6. 



III. CONCLUSIONS 

We have shown that the H A rational system related to the non-crystallographic root sys- 
tem H4 is exactly solvable with the characteristic vector (1, 5, 8, 12). This work complements 
the previous studies of the rational (and trigonometric) models, related with crystallographic 

nnn n n 

root systems (e.g. [6J - [9(, [3|]) and non-crystallographic root systems h(k) [10] and H 3 [11]. 



A certain significance of exploration of the H4 rational system is due to a fact that this model 
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is defined in four-dimensional Euclidian space. There are very few known exactly- solvable 
systems in this space - five-body Calogero-Sutherland (A4) and BC4 rational-trigonometric 
models among them. All of them are completely-integrable. 

Taking Coxeter invariants of H4 as coordinates provided us a way to reduce the rational iJ 4 
Hamiltonian to algebraic form. It gave us a chance to find the eigenf unctions of the rational 
i?4 Hamiltonian which are proportional to polynomials in these invariant coordinates. It 
seems correct that these eigenf unctions exhaust all eigenfunctions in the Hilbert space. It is 
worth noting that the matrix Aij(r) which appears in front of the second derivatives after 
changing variables in Laplacian from Cartesian to the H4 Coxeter invariant coordinates (see 
Eqs. (1121) ) has polynomial entries corresponding to flat space metric, hence the Riemann 
tensor vanishes. 

It should be stressed that it was stated in Lax pair formalism that the Hamiltonian 
of the rational system ([ID is completely integrable 12j. This implies the existence of 
three mutually-commuted operators (the 'higher Hamiltonians') which commute with the 
Hamiltonian forming a commutative algebra. It is known (see 1|) for the crystallographic 
systems that these higher Hamiltonians are the differential operators of the degrees which 
coincide to the minimal degrees of the root space (the Lie algebra) or their doubles for the 
An case. It may suggest that for the H4 rational system the commuting integrals might 
be differential operators of the orders 12, 20 and 30. Their explicit forms are not known 
so far. It seems evident that these commuting operators should take on an algebraic form 
after a gauge rotation (with the ground state function as a gauge factor), and a change of 
variables from Cartesian coordinates to the Coxeter invariant variables r's. Interesting open 
question is about a flag of invariant subspaces: would it be one of these two flags preserved 
by hn 4 ^ Following the experience with different integrable systems, it seems the integral(s) 
related with separation of variables do not enter to the commutative algebra. Therefore, the 
integral J 7 lays out of the commutative algebra of integrals. It might serve as an indication 
to a superintegrability of the H4 rational system. 

It should be pointed out that unlike the rational models of the crystallographic root 
spaces it is not possible to construct integrable (and exactly-solvable) trigonometric systems 
related to the non-crystallographic root spaces as a natural generalization of the Hamiltonian 
Reduction Method [1]. 

The existence of algebraic form of the H4 rational Olshanetsky-Perelomov Hamiltonian 
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makes possible the study of their polynomial perturbations which are invariant wrt the H± 
Coxeter group by purely algebraic means: one can develop a perturbation theory in which 
all corrections are found by linear algebra methods 13|. In particular, it gives a chance to 
calculate the if 4 Coxeter- invariant, polynomial correlation functions by algebraic means. 

Another important property of the existence of algebraic form of the rational Hamil- 
tonian is a chance to perform a canonical, Lie-algebraic discretization to uniform ([l^jl) 



and exponential |15| lattices, or mixed uniform-exponential lattices. In the case of all three 
lattices such a discretization preserves a property of integrability, polynomiality of the eigen- 
functions remains and it is isospectral. Making the weighted projective transformation (|18p 
of the H4 algebraic form (II ip we arrive at different algebraic form of the H4 Hamiltonian. 
Making then the Lie-algebraic discretization we arrive at a discrete model related to an 
original discrete model via change of variables. It can be considered as a definition of a 
polynomial change of variables for discrete operators. 



One can find the s/(2)-quasi-exactly-solvable generalization [16| of the H4 model which 
remains integrable. This is one of the first examples of quasi-exact-solvability related to 
non-crystallographic root systems. It complements the results obtained previously for all 



rational models related to crystallographic systems (see ll7|), for the hik) rational model 



and for H3 - each of these models admit a certain s/(2)-QES generalization in a 
form of sixth degree polynomial potential. 

Owing to the explicit knowledge of the ground state function ([2]) supersymmetric 



model can be constructed following a procedure realized in 18| for An rational model, in 



[191 ] for the BCn rational model and in 20] for the hik) rational model. It can be done 
elsewhere. 
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